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The averaging problem in cosmology is of fundamental importance. When applied to study 
cosmological evolution, the theory of macroscopic gravity (MG) can be regarded as a long-distance 
modification of general relativity. In the MG approach to the averaging problem in cosmology, the 
Einstein field equations on cosmological scales are modified by appropriate gravitational correlation 
terms. We study the averaging problem within the class of spherically symmetric cosmological 
models. That is, we shall take the microscopic equations and effect the averaging procedure to 
determine the precise form of the correlation tensor in this case. In particular, by working in 
volume preserving coordinates, we calculate the form of the correlation tensor under some reasonable 
assumptions on the form for the inhomogeneous gravitational field and matter distribution. We find 
that the correlation tensor in a Friedmann-Lemaitre-Robertson- Walker (FLRW) background must 
be of the form of a spatial curvature. Inhomogeneities and spatial averaging, through this spatial 
curvature correction term, can have a very significant dynamical effect on the dynamics of the 
Universe and cosmological observations; in particular, we discuss whether spatial averaging might 
lead to a more conservative explanation of the observed acceleration of the Universe (without the 
introduction of exotic dark matter fields) . We also find that the correlation tensor for a non-FLRW 
background can be interpreted as the sum of a spatial curvature and an anisotropic fluid. This may 
lead to interesting effects of averaging on astrophysical scales. We also discuss the results of averaging 
an inhomogeneous Lemaitre-Tolman-Bondi solution as well as calculations of linear perturbations 
(that is, the backreaction) in an FLRW background, which support the main conclusions of the 
analysis. 

I. INTRODUCTION 

The Universe is not isotropic or spatially homogeneous on local scales. The correct governing equations on cosmo- 
logical scales are obtained by averaging the Einstein equations of general relativity (GR). An averaging of inhomoge- 
neous spacetimcs can lead to dynamical behavior different from the spatially homogeneous and isotropic Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model in particular, the expansion rate may be significantly affected 0. 
Consequently, a solution of the averaging problem is of considerable importance for the correct interpretation of 
cosmological data. Unfortunately, this is a very difficult problem. 

There are a number of theoretical approaches to the averaging problem 0,111,1,0,0. In the approach of 
Buchert @, a 3+1 cosmological spacetime splitting that depends on the spacetime foliation is utilized and, in addition, 
only scalar quantities are averaged (and hence, in general, the equations are not closed and consequently additional 
assumptions are necessary). The perturbative approach [2, 8] involves averaging the perturbed Einstein equations; 
however, a perturbation analysis cannot provide detailed information about the averaged geometry. 

In all of these approaches an averaging of the Einstein equations is performed to obtain the averaged field equations. 
The macroscopic gravity (MG) approach to the averaging problem in cosmology is an attempt to Q give an exact (and 
tensorial) prescription for the correlation functions which inevitably emerges in an averaging of the field equations 
(without which the averaging simply amounts to definitions of the new averaged terms) . In MG the Einstein equations 
on cosmological scales with a continuous distribution of cosmological matter are modified by appropriate gravitational 
correlation (correction) terms. The MG approach provides a covariant method of averaging tensors, consequently it 
does not rely on assumptions regarding the nature of perturbations (e.g., in principle there are no approximations 
and no higher order terms are dropped). We shall adopt the MG averaging approach in this paper. 

The spacetime averaging procedure adopted in MG is based on the concept of Lie-dragging of averaging regions 42] , 
which makes it valid for any differentiable manifold with a volume n-form, and it has been proven to exist on an 
arbitrary Riemannian spacetime with well-defined local averaged properties [3j . Averaging of the structure equations 
for the geometry of GR leads to the structure equations for the averaged (macroscopic) geometry and the definitions 
and the properties of the correlation tensor. The averaged Einstein equations for the macroscopic metric tensor 
together with a set of algebraic and differential equations for the correlation tensors become a coupled system of the 
macroscopic field equations for the unknown macroscopic metric, correlation tensor, and other objects of the theory. 
The averaged Einstein equations can always be written in the form of the Einstein equations for the macroscopic 
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metric tensor when the correlation terms are moved to the right-hand side of the averaged Einstein equations, and 
consequently can be regarded as a geometric modification to the averaged (macroscopic) matter energy-momentum 
tensor 

MG is a non-perturbative geometric field theory with a built-in scale related to the spatial scale over which averages 
are taken (we recall that, in principle, averaging is performed over a 4- volume region). The microscopic field to be 
averaged is supposed to have two essentially different variation scales, A and Lh, satisfying A << Lh, where Lh is 
the horizon size related to the inverse Hubble scale. An averaging region must be taken of an intermediate size L 
such as A << L « Lh, so that the averaging effectively smooths out all the variations of the microscopic field of 
the scale A. In cosmological applications A is taken to be the scale on which astrophysical objects such as galaxies or 
clusters of galaxies have structure, and the size of the averaging space regions has been tacitly assumed to be ~ 100 
Mpc, or a fraction of the order of the inverse Hubble scale, and thus any terms (e.g., a cosmological constant or a 
curvature term) appearing in the correlation tensor might be expected to be related to the inverse Hubble scale. In 
principle the scale, given by the size of the spacetime averaging region, is a free parameter of the theory. 

A procedure for solving the MG equations with one connection correlation tensor was discussed in Q . The macro- 
scopic field equations were written in the form of the Einstein equations of GR, with a 'modified' stress-energy tensor 
consisting of the averaged microscopic stress-energy tensor l^( mlcro )^ anc j an additional effective stress-energy tensor 
C arising from the correlation tensor Z ||. In Q it was found that the averaged Einstein equations for a spatially 
homogeneous, isotropic macroscopic spacetime geometry has the form of the Einstein equations of GR for a spatially 
homogeneous, isotropic spacetime geometry with an additional spatial curvature term (i.e., the correlation tensor C 
is of the form of a spatial curvature term) . 

Therefore, assuming spatial homogeneity and isotropy on largest scales, then the inhomogeneities affect the dynamics 
through correction terms (the correlation tensor) of the form of a curvature term Q , which will dominate at late times 
and on largest scales. Thus even for FLRW backgrounds it is important to understand how these correction terms 
affect cosmological observations. For example, a spatially averaged metric is not a local physical observable: the 
averaged value of the expansion will not be the same as the expansion rate of the averaged geometry, because of the 
non-linear nature of the expansion. 

The spacetime averaging in MG utilizes bilocal averaging operators. The MG averaging scheme is especially simple 
in a proper coordinate system in which the bilocal operators takes on the simplest possible forms. In particular, 
any proper coordinate system is necessarily a volume-preserving (system of) coordinates (VPC), and in a pseudo- 
Riemannian spacetime the spacetime averages defined in proper coordinates are Lorentz tensors exactly like the 
averages in Minkowski spacetime; that is, VPC on an arbitrary differentiable metric manifold is a natural counterpart 
of the Cartesian coordinate system on a Minkowski manifold. A brief review of the spacetime averaging scheme 
adopted in macroscopic gravity and the role of proper systems of coordinates is presented in Appendix A. 

Spherically symmetric cosmological models are of special cosmological importance, partially motivated by the 
observed isotropy of the CMB. Therefore, it is important to study the averaging problem comprehensively within the 
class of spherically symmetric cosmological models (i.e., to determine the form of the MG equations in the case of 
spherical symmetry). We shall take the microscopic equations and effect the averaging procedure to determine the 
precise form of the correlation tensor in this case. 

In the next section we shall calculate the form of the MG equations in the case of spherical symmetry. The first 
step is to choose an appropriate spherically symmetric VPC system. It is also instructive to investigate the FLRW 
metric in VPC. We then make some reasonable assumptions on the form of the inhomogeneous gravitational field and 
matter distribution, and in section III we calculate the resulting form of the correlation tensor in both a FLRW and 
non-FLRW background. 

In section IV we average an inhomogeneous Lemaitre-Tolman-Bondi (LTB) solution. The first step in this calculation 
is to rewrite the LTB dust solution in VPC (which is done in Appendix B). This solution is presented as a perturbation 
about a flat FLRW model. In section V we then assume a spatially homogeneous and isotropic background and discuss 
the effect of perturbations (that is, the backreaction) on this FLRW background. 

In section VI we discuss the results obtained in light of recent observations, with particular emphasis on the effect 
of inhomogeneities on the local expansion rate. We discuss whether inhomogeneities and spatial averaging might lead 
to a more conservative explanation of the observed acceleration of the Universe (without the introduction of exotic 
dark matter fields). The conclusions are given in the final section. In Appendix C we briefly discuss the relationship 
between our work and the work of Buchert [6[ in the case of spherical symmetry. 

II. SPHERICAL SYMMETRY 

We shall calculate the form of the MG equations in the case of spherical symmetry. That is, we shall take the 
microscopic equations and effect the averaging procedure to determine the precise form of the correlation tensor, C a b , 
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in this case. 

We begin by choosing an appropriate coordinate system. Starting from the general form of the spherically symmetric 
metric, we first choose a new angular coordinate, u — cos(8), to eliminate any angular dependence in y/—g, where 
g = det(g a b). Next, we use the remaining coordinate freedom to set y/—g = 1; this is done by choosing an appropriate 
form for the 'radial' metric function that multiplies the spherical line element ds 2 (u, 4>). The line element is thus [43| 



ds 1 



-Bdt 2 



Adr 1 
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where the functions A and B depend on t and r. These are volume preserving coordinates (VPC) for the spherically 
symmetric metric [4|. It is the adoption of VPC that enables us to calculate the averaged quantities in a relatively 
straightforward manner. 

We next calculate [1(3] the form of the Einstein tensor G a b (note the position of indices) 
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and G t r = —AB~ 1 G r t , G , = G u u . The subscripts r and t on the metric functions A and -B denote partial differenti- 
ation with respect to r and t, respectively. 

We note that all terms in the expressions for G r t and G u u originate from the basic metric functions g rr and gu 
(A and B). The terms of the form V AB in G t t and G r r arise from derived 'radial' metric functions (e.g., the term 
V AB in equation ([2]) arises as a product of the metric components g u , g rr , guu and g^; also see the comment in 
subsection III. A. 3). In order to avoid unnecessary complications in the averaging procedure we shall not deal with 
derived quantities directly [44| . Therefore, we eliminate the s/AB terms by considering 
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and use the contracted Bianchi identity to determine the remaining component of the Einstein tensor. Taking averages, 
we now obtain the appropriate form for the MG equations and hence the correlation tensor C a b . For example, we 
have that 14511 
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In the above, by virtue of VPC, the average is now a simple average defined by 



(f(r,t)) = — 



d£ 



dr'f{r + r',t + t'), 
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which, for smooth functions with a sufficiently slowly varying dependence on cosmological time, essentially reduces 
to a spatial average in terms of the averaging scale L (see the next section) . 



4 



It is instructive to consider the FLRW metric in VPC. The metric is given by (JTJ) , with 



A-£- B-±- 



(9) 



where R = R(t) and F — F(r), subject to 4^ = (yl — kF 2 )F 2 , and k = —1,0 or 1. The Einstein tensor in these 
coordinates is given by 
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The spatial curvature term is given by G a b = — kRg diag[3, 1, 1, 1] (whereby spatial curvature term we mean the 
Einstein tensor corresponding to a spacetime with constant spatial curvature), with effective equation of state p e g- = 



III. INHOMOGENEOUS SPACETIMES 

The form of the correlation tensor now depends on the assumed form for the inhomogeneous gravitational field and 
matter distribution. Let us assume that 
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B(r,t) = (B( r> t)) 
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where r is a radial variable (and, strictly speaking, we are assuming that r < L). The assumptions (jlip and (|12|) 
constitute a spatial Fourier decomposition of the metric functions in which the variation in the timclikc direction is 
assumed small and the dominant source of inhomogeneity arises from a spatial variation of the gravitational field 
(thus the 4- volume average effectively reduces, in this case, to a smoothing on a spatial domain). Note that the 
coordinates t and r appearing in |T]) are not the usual 'time' and 'radial' coordinates; however, the unit magnitude 
timelike coordinate basis vector has zero vorticity, which implies the existence of a foliation of spacetime (where the 
r coordinate parameterizes the spatial hypersurfaces). Since the coordinate basis vectors dt and d r are independent 
(i.e., the metric is diagonal), it follows that variation along timelike and spatial directions is not coupled. Although 
other forms for the inhomogeneous gravitational field are possible (i.e., different assumptions to (|llll2p ). it is not 
expected that the main conclusions obtained in this paper will be qualitatively affected. In equations (jllj) and \12\ , 
L = Lq is treated as a parameter and in the calculations that follow Lq is effectively taken to be a small dimensionless 
parameter after a renormalization of the variables using the speed of light (set to unity) and the present value of the 
Hubble parameter, Hq, and a redefinition of the functions a n , b n . 

With these assumptions [46| we have essentially assumed that averages effectively become space averages with 
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We note that integrating the left-hand side of the first of these equations (i.e. , taking spatial averages) yields (A(r, t)) . 
We also note that 
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and 
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Thus the assumed forms for the inhomogeneous functions A and B satisfy a set of appropriate and self-consistent 
conditions. 

We shall expand in powers of L = Lq < 1. Since the coefficients (01,02, for example) are evolving functions of 
time, these expansions may only be valid for a transient period of time. We calculate the correlation tensor C a b up to 
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The 0(L 2 ) terms have been calculated, but we have not explicitly displayed them here. 



A. Lowest Order Calculation 



From eqns. (7) and (fl4|) we obtain, for example, C r t = + O(L), and from eqns. (5) and (6) 
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where C = C r r and 



£(t) = — [(01 - 3&0(oi + hi) + (ox - 36i)(oi + 61)] . 

We calculate C from the contracted Bianchi identities. We note that if C a & is isotropic (i.e., of the form of a perfect 
fluid) then C = ™ and C a & is of the form of a spatial curvature term. 

1. Bianchi Identities 

For the metric ([I]) and correlation tensor ([20]) . to 0(L°) the contracted Bianchi identities yield: 
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where C r = ^ and C' t = The solution of these equations (for C) depends on whether B r is zero or not. 

2. FLRW background 

In the case that B r = 0, as in the case of a FLRW background, equation (f2Tj) immediately yields C = and 
(A) tr = 0, and equation (|2"2")l then yields -^y = InBr 2 (where £q is a constant that depends on the spatial averaging 
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scale). Therefore, in this case we obtain 



C a b — £qR 



3 

10 

10 

1 



and C a b is necessarily of the form of a spatial curvature term. 

In order to fully reconcile (|2T))) and (|2"3"|) , we note that since A — R 2 F~ A in the FLRW case, we obtain 
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Here (A) must be interpreted as an averaged spatial curvature, so that (F 4 (r)) is replaced by (2F 2 F 2 + F 3 F rr ), which 
is constant. That is, the term (A) in the above equation must be interpreted correctly. 



3. Non-FLRW background 
If B r =^ 0, then eqn. (|2"Tj) can be integrated to obtain 
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(we note that, in general, this expression is different to what we would have obtained if we had averaged eqns. |T]), 
(TlTj) and (12"]). with the y/AB term, directly). Eqns and then yield 
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Since A = A(r,t) (B r ^ 0), in general the solution of this equation yields £ = £q (constant) and f(t) = g(r)A 2e , so 
that if / 7^ 0, A(r,t) is separable. In the latter case, in general B(r,i) is also separable, and the two separate terms 
in ([25|) can be of a comparable form. 

We note that C a b is necessarily anisotropic (and cannot be formally equivalent to a perfect fluid). We can always 
write 
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where II = — | <?(?") (AS) 1 / 2 — j- The correlation tensor C a b automatically satisfies the contracted Bianchi iden- 
tities (|2"Tj) and (f2"2"|) . It can be interpreted as the sum of a perfect fluid and an anisotropic fluid (when B r ^ 0). For an 
anisotropic fluid in spherically symmetric coordinates the energy-momentum tensor is of the form diag [— /i, p\ \ , p± , px] , 
where p\ | = p+ ^ir and p± = p~ ^tt, and ir is the anisotropic pressure. From above, we see that if the (total) correlation 
tensor C a b is interpreted as an anisotropic fluid, it follows that II = — ir and /i = 3p. 

If both terms separately satisfy the contracted Bianchi identity, then the first term can be interpreted as a spatial 
curvature term. The second term can be interpreted as an anisotropic fluid with p± = and p|| = — p„a (which is 

similar to the equation of state for a cosmological constant). For f(t) = g{r)A (g(r) ^ 0) and A = F {r)R (t), 
we obtain pug = 2£qF (t)R~ 2 = g{r)F ' R^B 1 / 2 , so that B is separable and of the form B = b{r)R (t) (compare 
with eqn. (j9|). One solution gives p e Q = 0, so that in this case C a & is of the form of a perfect fluid and is thus 
necessarily of the form of a spatial curvature term. 

Although the correlation tensor C a b satisfies the contracted Bianchi identities, when interpreted as the sum of a 
spatial curvature perfect fluid and an anisotropic fluid through (|26|) . the two separate fluid do not in general satisfy 
separate conservation equations. However, eqn. (|22p can be rewritten in the form of a conservation law for the 
anisotropic pressure n, 
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in VPC where the metric is given by eqn. (Q} and u a is comoving to order 0(L 2 ) (note that the expansion and shear 
are given by 

1 * '3-i + -± , (28) 



2Bi \ B J ' 2x/6Bi V ^ B 
respectively; compare this with eqn. (|72| in Appendix C). 



^. Anisotropic Fluid 

The second term in the correlation tensor is of the form of an energy-momentum tensor for an aniso trop ic fluid, 
with energy density fi, a pressure p|| parallel to the radial unit normal and a perpendicular pressure v 1 [47] . Fluids 
with an anisotropic pressure have been studied in the cosmological context for a number of reasons: an energy- 
momentum tensor of this form arises formally if the source consists of two perfect fluids with distinct four- velocities, 
a heat conducting viscous fluid under some circumstances, a perfect fluid and a magnetic field, and in the presence of 
particle production [ill [T2I [l3j. In particular, the energy-momentum tensor of a cosmic string [lij is of the form of 
an anisotropic fluid with /i = —p\\, p± = (such an equation of state also arises in other early universe applications). 
Anisotropic fluids in spherically symmetric cosmological models have been studied in pH fl2l. [l5| . In addition, the 
energy-momentum arising from the gravitational field of a global monopole is formally an anisotropic fluid which is 
static and spherically symmetric [Tij . We also note that in an investigation of the consequences of an imperfect dark 
energy component on the large scale structure, the effect of anisotropic perturbations (due to the dark energy) on the 
cosmic microwave background radiation was studied. It was found that an anisotropic stress is not excluded by the 
present day cosmological observations [l6l |. 

Let us comment on the astrophysical applications of an anisotropic fluid. It is known that dark matter is a major 
constituent of the halos of galaxies fl7T ] . By an analysis of observed rotation curves, under reasonable assumptions 
(e.g., that galaxies can be modeled as spherically symmetric) it has been found that the dark matter is of the form of 
an anisotropic fluid [l|[ . This has been taken up in , in which the consequences of anisotropic dark matter stresses 
are discussed in the weak field gravitational lensing (where it was noted that in any attempt to model dark matter in 
galactic halos with classical fields will lead to anisotropic stresses comparable in magnitude with the energy density) . 

Finally, we note that for a 4-dimensional spacetime with a metric of the form g a b — diag[— 1, 1, (1 + cu 2 ) -1 , (1 + cu 2 )] 
(in VPC), where c is a constant, we have that G a & = diag[c, c, 0, 0], The metric is of the product form R 2 x S 2 , 
and is therefore the tensor product of a 2-dimensional flat space and (for c < 0) a 2-dimensional sphere (which are 
two 2-dimensional spaces of constant curvature). Hence, the second term can also be interpreted in terms of spatial 
curvature (although we again note that each of the two terms, namely the spatial curvature term and the anisotropic 
term, do not separately satisfy the contracted Bianchi identity). 



B. Further calculations 

1. Higher order terms 

We can consider the contribution of the higher order 0{L) and 0(L 2 ) terms. In the case of a FLRW background 
(with B r — 0, -A\ = IqR~ 2 , where £ is constant), the O(L) terms all vanish. We can see from Appendix B that this 

occurs trivially when (a\b\ — aib\) = 0, which occurs for a sine series expansion (only) or a cosine series (only), or for 
a single trigonometric series with a composite argument (using the double angle formulae). The 0(L 2 ) contributions 
to C r t then vanish, and the 0(L 2 ) contributions to the remaining components of the correlation tensor are of the 
form of a spatial curvature term (the 0{L 2 ) terms are not displayed explicitly in the Appendix). That is, the effect 
of the higher order terms is simply to renormalize the spatial curvature term. 

In the non-FLRW background case, in general we must have (ai&i — aib\) = 0, and the higher order terms do not 
play any significant role (as above). 



2. Discussion 



In equations ((11)) and (fT2|) . L is essentially treated as a dimensionless parameter, which is sufficient to the lowest 
order of approximation. In principle, in the cosmological setting L depends on the Hubble scale H~ 1 and might also 



8 



be related to a scale dependent on structure formation, both of which vary with cosmological time. Therefore, in 
general, L will be time dependent. 

To lowest order we assume that L = Lq is fixed and integration is taken over a comoving domain (and presently 
Lq ~ 1CP 1 ). Assuming L is time dependent, L = L(t), we have that a typical correction term is of the form 



1 1 



T/2 



Assuming time evolution is of the order of the Hubble scale, we have that Ct ~ 0(Lq x HqT), where Hq is the 
current value of the Hubble parameter. Clearly such corrections are of order O(L ) compared to the contributions 
calculated above. Moreover, these corrections are negligible over small time averaging scales T (compared to Ho; i.e., 
H T small). In addition, we have that 

where 

C d T=(^{f-(f)}}~L HoU, (31) 

where the term C in the particular case of inhomogcncities of the form (TTT j) / (fT2 |) is negligible. Therefore, TCqt ~ 
O(Lq) x HqT. These corrections are consequently of higher order and generally will only renormalize the spatial 
curvature term. 



3. Other inhomogeneous models 

The form of the correlation tensor depends on the assumed form for the inhomogeneous gravitational field and 
matter distribution. We could consider alternatives to the form of the inhomogeneous metric (fTT j) / (fT2 |) . We shall 
consider two alternative approaches here. First, we shall average an exact inhomogeneous Lemaitre-Tolman-Bondi 
solution. Second, we shall discuss a linear inhomogeneous perturbation of an exact FLRW model. 

However, the main conclusions of this section will not be affected; namely, in most applications of interest the 
correlation tensor is of the form of a spatial curvature, but in general it is not even of the form of a perfect fluid. 
Moreover, higher order corrections are not expected to lead to significant effects; e.g., they alone cannot account for 
a current acceleration. 



IV. lemaItre-tolman-bondi model 



Let us consider av eraging an exact solution. The spherically symmetric dust solution is the exact Lemaitre-Tolman- 
Bondi (LTB) model [ll|,|20j], which can be regarded as an exact inhomogeneous generalization of the FLRW solution. 
In the dust LTB model, from the Gauss-Codazzi equations the Einstein tensor has the form of a spatial curvature 
tensor on spacelike hypersurfaces (which we recall is not the same as the projected Einstein tensor). Various aspects 
of the averaging problem in LTB spacetimes have been studied [2l|, |22| . 

The first step is to take the LTB solution [HI, [HI and rewrite it in VPC. This is done in Appendix B. We note that 
taking averages using VPC is of interest in its own ri ght, and is an advantage in that averaging can now be done in 
both space and time (this will be discussed further in [23j ). From the Appendix, we obtain 



2- 4 dtd X + ^ 



R 2 



du 2 



+ (1 - u 2 )d(f> 2 



(32) 



The constraints of the original LTB metric become a defining equation for U(t,x) (|54|) and a differential equation for 
R(t, x) (|55|) . which then ensures an exact dust solution with density G u (t, x). Since in VPC the velocity of the dust 
flow is u a = (l,U(t,x),0,0), the Einstein tensor components satisfy G tx = UG U ,G XX = U 2 G U ,G UU = G^ = (see 
eqns. (|5f?)) ). 

In appendix B we explicitly construct the FLRW dust models in VPC. The spatially flat (Eq = 0) FLRW model 
in VPC is given by (|57|) . The spatially closed (Eo < 0) FLRW model in VPC is given in eqn. (|6Tj) in terms of an 
expansion (of R, U and G M ) about the spatially flat FLRW model with E = (a similar expression exists for the 
E Q >0 FLRW model). 
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A. A Perturbative Solution 



We shall assume that tB{f) is zero, which implies that the bang time is uniform and we are consequently restricting 
our choice of LTB models to those with no decaying modes. Such models are of interest at later times, and particularly 



in the study of structure formation [24(, and are suitable for our purposes here. We shall also consider solutions of 
the LTB metric in VPC as perturbations about the spatially flat FLRW model given in {57}. In this respect our 
approximate solution will be an expansion with respect to Eq and we require the Einstein tensor to have the form of 
([56]) (i.e., the form of dust, after truncation of terms of 0(Eq) or higher). We begin by making the following ansatz 
on the form of R 

R{t, x) = Ro + a lX a t b E + a 2 x c t d El, (33) 

where a.\, «2, a, 6, c and d are constants to be determined from requiring the Einstein tensor has the form of a dust 
solution up to order Eq. Substituting (|54p into ([55| gives a partial differential equation (PDE) involving only R, 
a subsequent substitution by (f3"3"|) then shows that the first non-trivial term in the PDE occurs at O(Eq). At this 
stage we are only interested in a perturbative solution of the PDE, therefore to obtain necessary conditions for a dust 
solution we simply require that the coefficient of Eq vanish, which yields three cases (a = 1/3, —2/3 or 1/3 — b/2). 
Choosing a, we can use eqns. (|33|) and ([54]) to obtain U(t,x). Calculating the Einstein tensor and requiring it have 
the form of (fS"6| allows us to determine the remaining constants (the details are discussed in [231]). 

If a — 1/3, a number of subcases occur. A typical solution (e.g., b = 0, c = 5/6 and d = —1) gives rise to 
G tt = 4/(3i 2 ) - 2aia;- 2 / 3 £ , o + 0(E$) and G^ = 0(E$), so that the truncation of Eq and higher terms results in 
an Einstein tensor of the form of dust ([56]) . Other solutions give rise to a G tt with no Eg terms (but containing 
higher orders of Eq, with G^ beginning at these higher orders) and G u and G xx components with no 0(1) terms 
(and beginning with Eq terms, whereas the other components begin at higher orders of Eq). After averaging the 
(— 2q!ix~ 2 / 3 Eq) contribution, we obtain a correction term to the density which is independent of t. This constant 
correction in the dust model may be related to a cosmological constant or to an anisotropic source (this will be further 
investigated in [23[). 

If a = —2/3, then the resulting dust solutions are of the form of a flat or curved FLRW model (the closed model is 
given by (|59p). That is, the solution in this case is typically of the form of a flat FLRW model with spatial curvature 
corrections. 

If a = 1/3 — b/2 (where b ^ 0,2, these cases are discussed above), setting c = —1/3 — 6, d = 2(6 + 1) and 
a 2 = — 81a 2 6(6 — 2) 2 /(80(36 — 1)) gives a dust solution up to order Eq. Unlike previous cases, here we have an 
arbitrary power of x, and there are two free parameters b and a,\] however, these solutions are not more general 
because other free parameters are constrained. 

We can obtain more general perturbative dust solutions containing more free parameters by a superposition of the 
solutions discussed above. For example, one such solution is found by adding a solution with a = 1/3, 6 = 0, c = 5/6 
and d = —1 and a solution with 6 = —16/3, a = 3, c = 5, d = —26/3 and a 2 = — 1452a 2 /85, to obtain a solution of 
the form (to order Eq) 

R(t, x) = a R Q + [aix 1 / 3 + ^xH-^'^Eq + [a^^r 1 - (1452/85)/3 2 x 5 ^ 26 / 3 ]£ 2 , (34) 

with corresponding forms for U(t,x) and G'*, thus giving rise to another (more general) perturbative dust solution. 
Therefore, we can construct perturbative LTB solutions which can be interpreted as having both spatial curvature 
and constant correction terms. 



V. COSMOLOGICAL PERTURBATIONS 



A. Backreaction 



The theoretical approach is to solve the full problem to obtain the equations satisfied by the averaged quantities, 
without assuming a given background. An alternative but more practical approach is to assume a spatially homo- 
geneous and isotropic background and study the effect of perturbations (that is, the backreaction) on this FLRW 
background [3, H, [25|, [26[ . The starting point is the Einstein equations in an appropriately defined background [25l \2l\ . 
The Einstein and energy-momentum tensors are then expanded in metric and matter perturbations up to second 
order. The linear equations are assumed to be satisfied, and the spatially averaged remnants provide the new back- 
ground metric which takes into account the backreaction effect of linear fluctuations computed up to quadratic order. 
The backreaction has been studied for scalar gravitational perturbations [25| . and it was found that the equation of 
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state of the dominant infrared contribution to the energy-momentum tensor which describes backreaction can take 
the form of a negative cosmological constant. This has led to the speculation that gravitational backreaction may 
lead to a dynamical cancelation mechanism for a bare cosmological constant [28l |. Since, in the perturbative approach, 
the averaged Einstein tensor will rapidly come to dominate over the correlation tensor, it has been argued that this 
might also explain the presence of a source of late-time acceleration [28| . 

The main aim is to investigate the effect of these perturbations on the local expansion rate and to see how it might 
differ from the background expansion rate. In an important study, using a 3+1 split and the Zeldovich approach 
and assuming inhomogeneous perturbations about a dust Einstein de Sitter (FLRW) background, the equations 
governing the time dependence of the scale factor due to backreaction were obtained by spatial averaging 3. The 
metric perturbations were assumed small, even when the density contrast is large (much larger than unity). It was 
found that the scale factor dependence on the correlation terms acts like a (negative) spatial curvature term (and, 
curiously, that the age is greater than in the exact flat background FLRW model) Q. Typically perturbations are 
small corrections, but since they are time dependent they can become larger, although likely vanishing asymptotically 
to the future [2f3]. 

Therefore, the resulting correlation terms are simply of the form of a spatial curvature term in the linear perturbation 
analysis. This is true in general, and is certainly true in spherically symmetric cosmological models. There are, as 
mentioned earlier, problems with the perturbative approach. First, the perturbation scheme breaks down when 
perturbations become significant and affect the background. Second, there are potentially gauge effects arising from 
the choice of hypersurface on which to do spatial averaging. 

B. Discussion 

Recent observations are usually interpreted as implying that the Universe is very nearly flat, currently accelerating 
[30l | and indicating the existence of dark matter and dark energy [29| . A cosmological constant or a negative pressure 
fluid (or quintessence field) are candidates for the dark energy. However, as noted earlier, inhomogeneities can affect 
the dynamics and may significantly affect the expansion rate of the spatially averaged "background" FLRW universe 
(the effect depending on the scale of the initial inhomogeneity) 0]. Therefore, a more conservative approach to 
explain the acceleration of the Universe without introduction of exotic fields is to utilize a backreaction effect due to 
inhomogeneities of the Universe. 

It has been suggested that backreactions from inhomogeneities smaller than the Hubble scale could explain the 
apparently observed accelerated expansion of the Universe today. This has been investigated by studying the effec- 
tive Friedmann equation describing an inhomogeneous Universe after smoothing out of the sub-horizon cosmological 
perturbations, and it has been suggested that the acceleration in our Hubble volume might be possible even if local 
fluid elements do not individually undergo accelerated expansion (3ll [331] - However, in [34| it was claimed that the 
perturbative effect proposed amounts to a simple renormalisation of the spatial curvature, and in other work it has 
been argued that the acceleration cannot be explained by the effects of inhomogeneities [111 [36|. However, after 
density fluctuations in the Universe grow to be non-linear and begin to re-collapse, the perturbative expansion breaks 
down and reliable results cannot be obtained beyond this based on perturbative calculations. More recently, a solution 
using the gravitational backreaction of long wavelength (super-Hubble) fluctuation modes on the background metric 
was presented [28[ , and it was shown that in the presence of entropy fluctuations backreaction of the non-gradient 
terms is physically measurable (compare with [36[). 

In further work [26| the relationship between backreaction and spatial curvature using exact equations which do 
not rely on perturbation theory was studied in more detail, and it was argued that even though the effect does not 
simply reduce to spatial curvature, the acceleration that results is accompanied by a growth of spatial curvature to 
an extent that it is unlikely to be compatible with the CMB data. On the other hand, an explicit example of an 
inhomogeneous Universe has been presented that leads to accelerated expansion after taking spatial averaging [22| . 
The model universe is the LTB solution and contains both a region with positive spatial curvature and a region with 
negative spatial curvature. It was found that after the region with positive spatial curvature begins to re-collapse, 
the deceleration parameter of the spatially averaged universe becomes negative and the averaged universe starts 
accelerated expansion. Further examples, in which the assumption of spherical symmetry is relaxed, are discussed in 
[37l ]. In addition, inhomogeneities can lead to a reinterpretation of the luminosity distance of cosmological sources 
in terms of its redshift, which may account for the observed acceleration [38j]. However, it should be reiterated that 
there are subtleties when dealing with spatially averaged quantities, even if the spatial averaging is over a limited 
domain, and that the results discussed above may not apply to the quantities of physical interest [37[ • This point has 
also been further emphasized in [3(| (also see [28l l33l] ) . 
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VI. CONCLUSIONS 



We have calculated the form of the MG equations in the case of spherical symmetry. By working in VPC, we 
calculated the form of the correlation tensor under some reasonable assumptions on the form for the inhomogeneous 
gravitational field and matter distribution. The main result of this paper is that the correlation tensor in a FLRW 
background must be of the form of a spatial curvature, while for a non-FLRW background (with B r 7^ 0) the correlation 
tensor can be interpreted as the sum of a spatial curvature and an anisotropic fluid. We note that working in VPC (in 
which the determinant of the metric and hence the volume element is constant) is useful in its own right, particularly 
in the context of averaging. 

The cosmological result that in the spherically symmetric case the averaged Einstein equations in an FLRW back- 
ground has the form of the Einstein equations of GR for a spatially homogeneous, isotropic macroscopic spacetime 
geometry with an additional spatial curvature term, confirms the results in previous work in which we were able to 
explicitly solve the MG equations to find a correction term in the form of a spatial curvature @ . The results of the 
calculations regarding averaging of an exact inhomogeneous LTB solution (presented above), as well as calculations of 
linear perturbations (that is, the backreaction) in a spatially homogeneous and isotropic background and the results 
of Buchert [(| also confirm and support this result . 

The MG method adopted here is an exact approach in which inhomogeneities affect the dynamics on large scales 
through correction terms 48 1. Averaging can have a very significant dynamical effect on the evolution of the Universe; 
the correction terms change the interpretation of observations so that they need to be accounted for carefully to 
determine if the models may be consistent with an accelerating Universe. Averaging may or may not explain the 
observed acceleration. However, it is clear that it cannot be neglected, and a proper analysis will not be possible 
without a comprehensive understanding of the affects of averaging. 

On cosmological scales (of the order of the inverse Hubble scale), in which B r = and we have a FLRW background, 
averaging only gives rise to a spatial curvature term. However, the effects of averaging on astrophysical scales, such as 
galactic scales, are also of interest. Under the assumption that a galaxy can be approximated as spherically symmetric, 
where the background has B r 7^ 0, averaging is found to give rise to a correlation tensor of the form of an anisotropic 
fluid. This is of particular interest since, as noted earlier, dark matter in galactic halos is more accurately described 
by an anisotropic fluid [18| . 

In this paper we have also discussed averaging in an inhomogeneous LTB solution. Writing the LTB solution in 
volume preserving coordinates, we found a perturbative solution in which the correlation term can be interpreted 
as the sum of a spatial curvature term and a constant correction term. We also discussed linear inhomogeneous 
perturbations (that is, the backreaction) on an exact FLRW background. It was noted that the resulting correlation 
terms are simply of the form of a spatial curvature term. 

If the underlying microscopic spacetime has positive spatial curvature (as perhaps suggested by recent observations 
[2^ . l30l . HH), then we could obtain a cosmological model which is 'closed' on local scales, but as a result of the 
MG correlations behaves dynamically on macroscopically large scales as a flat model, which might have considerable 
physical implications. Indeed, cosmological models which act like an Einstein static model on the largest scales are 
possible even for models with zero or negative curvature on small scales; thus at late times (and on the largest scale) 
a sp atial curvature term will dominate the dynamics and the correlations might stabilize the Einstein static model 
[39^ (also see Section 3.3 of [37|L 

The MG analysis presented here is a self-consistent analysis. However, the MG approach is quite complicated and 
can be difficult to implement in practice. Therefore, it is of interest to compare our results to the work of Buchert [3^ |. 
In this latter approach a 3+1 split of the equations is effected (which introduces some gauge issues that presumably 
can be appropriately dealt with in the cosmological setting). More importantly, only scalar quantities appear in the 
averaged equations. This implies some sort of 'truncation' of the Einstein equations in order for the equations to 
reduce to scalar equations. As a result, in general the Buchert equations are not closed. In the approach taken here, 
the actual averages are constructed and therefore, in principle, the forms for the averaged quantities take on a specific 
form. Consequently, our approach is more restrictive in the sense that the system of governing equations is closed 
and so no further assumptions to close the system, which may or may not be physical, are necessary. However, it is 
anticipated that it is possible to derive the averaged scalar equations of Buchert [g] as some appropriate limit of MG. 
Whether any significant effects are neglected in the Buchert approach could then be determined. 

In Appendix C we show that in the spherically symmetric case the governing Einstein equations can reduce to 
equations in terms of scalar quantities in some circumstances. Therefore, we may be able to compare our results 
with the work of Buchert Q in this case. Certainly, a spatial curvature term appears in Buchert's scheme when 
averaging in a FLRW background, and it is expected that our analysis is consistent with the work of Buchert in a 
more fundamental sense. Unfortunately, the perfect fluid case is very complicated for comparisons, and so it is more 
sensible to first try to make the comparison in the dust case (i.e., the LTB model). Some brief comments are made 
in Appendix C and we shall pursue this further in [23[. Indeed, we note that for the models studied in section III, 
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if the correlation tensor is of the form of a perfect fluid then necessarily it is of the form of a spatial curvature (and 
this result is thus trivially consistent with the work of Buchert). On the other hand, if the correlation tensor is of the 
form of an anisotropic fluid, then a comparison with the work of Buchert is not possible. 

VII. APPENDIX A: MACROSCOPIC GRAVITY 

Let us review the spacetime averaging scheme adopted in macroscopic gravity (MG) [3j. It is based on the concept 
of Lie-dragging of averaging regions and is valid for any differcntiablc manifold. Choosing a compact region ScM 
in an n-dimensional differentiable metric manifold (M., g a p) with a volume n-form and a supporting point x S E to 
which the average value will be prescribed, the average value of a geometric object, Pp(x), x € M. , over a region E 
at the supporting point x € E is defined as 

n{x) = r±rj plfoaOV^dV = <P^) , (35) 
where Vs is the volume of the region E, 

= j J~gd n x , (36) 

the integration is carried out over all points x 1 € E, g' = det(g a p(x')) and the bold face object Pp(x,x') is a bilocal 
extension of the object Pp(x), 

p a p{x,x') = W^{x,x') P ^{x')W^{x\x) , (37) 

by means of the bilocal averaging operator Wpi(x,x') and its inverse Wp (x',x). The averaging scheme is covariant 
and linear by construction, and the averaged object pp has the same tensorial character as pp". As a result of the 
coincidence limit (\im x '^ x Wgi(x,x ) = Sp) and the idempotency condition, the average tensor Pp(x) takes the same 
value as the original tensor p%{x) when the integrating region £ tends to zero, which implies that the averaging 
procedure commutes with the operation of index contraction. 

In order to obtain the averaged fields of geometric objects on M. it is necessary to assign an averaging region Yj x to 
each point x of U CM, where the averaging integral is to be evaluated. To calculate directional, partial and covariant 
derivatives of the averaged fields, regions are related by Lie-dragging by means of a second bilocal operator, which 
can also be taken to be Wp (a/ , x) (which satisfies a divergence-free condition in order for Lie dragging of a region to 
be a volume preserving diffeomorphism) Q. The commutation relations simplify, and the differential constraint for 
the idempotent bilocal reduces to 

W[l 7] + Wg' a ,w£ = , (38) 

which has the general solution 

Wp*'(x',x) = f?(x')f- li p(x) (39) 

where f° t {x)d a — fi is any vector basis satisfying the commutation relations [fi, fj] = C^fk, with constant structure 
functions C^. In any n-dimensional differentiable metric manifold (M., gafj) with a volume n-form there always exist 
locally volume-preserving divergence- free operators Wp (x',x) of the form (|3"9")) @. 

A. Proper systems of coordinates 

We can consider the MG averaging scheme for a particular subclass of operators in which the averages and their 
properties are especially simple. Such a coordinate system is the analogue for MG of the Cartesian coordinates in 
Minkowski spacetime [4| . Let us hereby restrict the class of solutions of the equations (|3"5|) to the subclass satisfying 
[fi, fj] = 0; that is, Cy = 0. In this case the vector fields f" constitute a coordinate system and there always exist n 
functionally independent scalar functions (j? (x) such that the vector and corresponding dual 1-form bases are of the 
form 
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Thus, the bilocal operator Wg (x',x) becomes 



dx a 



W ^ X ^ = W0^' (41) 

Since they are functionally independent, the set of n functions </> z (x)can be taken as a system of local coordinates on 
the manifold A4, which will be called a proper coordinate system [J]. Therefore, in a proper coordinate system the 
bilocal operator (V, a;) takes the simplest possible form 

= Wg(x',x) lxa =t* = = 4> ( 42 ) 

where the bilocal Kronecker symbol Sp is defined as 8^ = 8f Sp. 

The definition of an average consequently takes on a particularly simple form when written using a proper coordinate 
system. The existence of volume- preserving bilocal operators of this form was proven in |4s] . Moreover, any proper 
coordinate system with a corresponding divergence-free bivector is necessarily a system of volume-preserving coordi- 
nates. In the case of a pseudo-Riemannian manifold the Christoffcl symbols, r^ Q , vanish and partial differentiation 
and averaging commute in VPC. 

Consequently, if the manifold (A4, g a p) is a pseudo-Riemannian spacetime, the spacetime averages defined in proper 
coordinates are Lorentz tensors, precisely like the averages in Minkowski spacetime. The average value of a tensor 
field pp(x), x £ £, over a compact space region S and a finite time interval At at a supporting point (t, x a ) G At x S 
is thus 

(p a p(t,x a )) £ = ^- s i [p%(t + t',x a + x a ')dt'd 3 x' , (43) 



At JS 



where Vs is the 3-volume of the region S, which is usually taken as a 3-sphere of radius R around the point x a at the 
instant of time t, Vs = J s d 3 x'. 

One issue of concern in the MG approach is the question of uniqueness; to what extent do spacetime averages 
depend on the choice of the bilocal operator. In the context of the present analysis, this raises the question of whether 
the results obtained in this paper could depend on the choice of VPC (Q]). It is clearly of interest to study this question, 
and we hope to return to this in future work. However, it is strongly anticipated that the main conclusions of this 
paper will not be qualitatively affected by the choice of VPC; namely, that in physical applications the correlation 
tensor is of the form of a spatial curvature (while in general the correlation tensor is of the form of an anisotropic 
fluid). 



VIII. APPENDIX B: LEMAITRE-TOLMAN-BONDI MODELS 



A spherically symmetric solution of the Einstein equation with dust field is given by the Lemaitre-Tolman-Bondi 
(LTB) solution [1JJ with metric 



ds 2 = -dr 2 + 4M^dr 2 + R 2 dQ 2 , 



(«y 



2E(r) 2M(r) 



(44) 
(45) 



where E(r) and M(r) are arbitrary functions of r. The solution of eqn. (|45p can be written parametrically by using 
the variable rj = J dr/R, 



R (V, r) = ^ [l - cos (vQ£W V ) 



-2E(r) 



y/-2E(r) 



y/-2E(r)ri 



By introducing the following variables 



R(r,r) 2E(r) 6M(r) 
a(r,r) = , k(r) = — , p (r) = 3— 



(46) 



(47) 
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the metric and the evolution equation for the scale factor a(r, r) become 



ds 2 



-dT 2 



_ k(r) 



Po( r ) 



(48) 
(49) 



Ean. (|49p is the same as the Friedmann equation with dust, and we can regard the LTB solution as a model of an 
inhomogeneous universe whose local behavior is equivalent to a FLRW universe with a spatial curvature k(r). 
As a specific case, Nambu and Tanimoto [221 ] assumed the following spatial distribution of spatial curvature: 



fc(r) = -!j[20(r-ro)-l] 



< r < L, < r < L 



(50) 



and assumed that po(r) = po =constant. For < r < tq, the solution is that of a spatially open FLRW universe and 
for ro < r < L, the solution is that of a spatially closed FLRW universe. 



A. Volume Preserving Coordinate System 



Starting with the LTB metric in the standard coordinate system (r, r, 9, </>) above that is aligned with the fluid flow, 
we obtain a volume preserving coordinate system (VPC), (t, x, u, 4>), by making the following coordinate transforma- 
tion 



t = T, X 



R(t, r) 2 R r 



■.dr, u = cos ( 



Defining U(t,x) := x T , and regarding R = R(t,x), the line-element becomes 

du 2 



ds' 



1 



El 

R 4 



,9 U , , dx 2 
dt 2 - 2—dtdx + — + R 



1-u' 



+ (1 - u 2 )d(/) 2 



R 4 R 4 

which has g = — 1 as desired. The constraints on the original LTB metric ensuring a dust solution with density 



now become 



and 



U(t,x) 



'R(t, r) 2 R r ' 



2RtR x + RRt 
' 2R 2 +RR XX 



2(3R x R t - RR X U X )U - 2RR x V t + (7R 2 X + 2RR XX )U 2 - Rf - 2RR tt + R 2 X R 4 = 1 



(51) 



(52) 



(53) 



(54) 



(55) 



Using eqn. (|5l|) . we can view equation (|55|) as a differential constraint for R(t,x) (where U(t,x) is then derived once 
R(t, x) is known). As a result of eqns. (fSTj) . (|54|) and (f55|) . the Einstein tensor has the following form 



G' 



G TT {t,x), G L 



XJG l \ G x 



U 2 G U . G u 



0. 



(56) 



These equations describe the corresponding LTB dust solution in volume preserving coordinates. We note that the fluid 
was comoving in the original coordinate system whereas in VPC the velocity of the dust flow is u a — (1, U(t, x), 0, 0), 
and hence is no longer comoving but is normalized, u a u a = — 1. 



B. Friedmann-Lemaitre-Robertson- Walker Cosmologies 



In order to compare with calculations we present the FLRW dust models in VPC. In the original coordinate 
system the FLRW models have E(r) = E r 2 , M(r) = M r 3 and a constant bang time is. Setting L = \E \ 3 / 2 /M 
throughout, the spatially flat (E = 0) FLRW model in VPC is given by 
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2x 

R(t, x) = (3a;) 1 / 3 , U(t, x) = ——, (57) 

t — IB 



giving an Einstein tensor component (the other components are determined from (|56|) ) 



G " = w^s?- (58) 



The exact spatially closed (E < 0) FLRW model in VPC is given by 

s 1 . . , . . 6\/2Lqx sin ri 

R{t,x) = — =— sm9 (1 — cos 77), U[t,x) = -- — , 

2y/2L (l-cosr/) 2 

32\/2~L 3 a; 

7? - sinry = 2V2L (t - t B ), 9 - sin cos 9 - - x -" !! 



(59) 



(1 — cos?7) 3 ' 

where r\ = rj(t) and 9 = 9{t, x). The resulting Einstein tensor component is 

G« = - 48L " w (60) 

(1 — COS Tjj 

It is of interest to consider the form of R, U and G u in the closed FLRW model written as a perturbation about the 
spatially flat FLRW model. Expanding about Eq = and defining Rq — (3a;) 1 / 3 , we obtain 

= ijl + ^ - 2./ V /3 (19 , 3 ^/3 + 126(1 - tB ) % \ 

V l5M 2/ '\t-t B )^ 1575M 4/3 (t-t B )8/3 7 

^(t,^) - T — — H "777^0 h +U{E ), (52) 

5Mq/ (i-t S ) 1/3 175Mq /3 

4 2 ■ fi 2 / 3 1 02 • fi 1 / 3 

G " = «i \ u ^ r E ° + 7 E l + °( E o)- (63) 

Kt-ts) 2 5M 2/3 (t - t B ) 4 /3 175A/ 4/3 (t - i B )2/3 

In the above perturbation scheme we have assumed that I?o (and x and £) are small. The spatially open (Eq > 0) 
FLRW model in VPC can also be displayed (it is similar to the closed case but with hyperbolic trigonometric functions 
replacing trigonometric functions and appropriate sign changes- see [Hj|). 



IX. APPENDIX C: SPHERICALLY SYMMETRIC MODELS IN THE 1+3 FORMALISM 

Using the Uggla and van Elst 1+3 formalism (40| for a perfect fluid energy-momentum tensor, in the case of spherical 
symmetry we have the evolution equations 

9 = -l9 2 + (e 1+ u-2a)(u)-^ (a+) 2 - I (/,+ 3p) + A (64) 
cr'+ = - 6 a + - (ei + it + a) (u) - *S+ (65) 
a = -\{9 + a+){a + u) (66) 

2 K = - 2 -{9 + a + ) 2 K (67) 

A = - (P+P)0 (68) 
and the Friedmann constraint 

= 1 -e* + \*R- 1 -a%-»-k (69) 

(and the spatial constraint equations), where 2 K(t,x) is the 2-curvature of the spheres, u is specified by choosing a 
temporal gauge, and p is specified by the fluid model. 
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There is only a single shear component, a = er + , so that in the spherically symmetric case the shear can be 
described by a single scalar. The curvature, in some cases, can also be described by a scalar. Therefore, in the 
spherically symmetric case we may be able to use scalar equations in some appropriate limit to describe the model. 
In particular, choosing a gauge in which u — and defining 

R d = 2 K 

Q d = -\« 2 

we obtain from the Friedmann constraint 

° = \° 2 + \ Rd + \ Qd ~ M ~ A + 2ei(a) ~ 6 ° 2 - 
The evolution eqns. (|65|) and (|67|) yield 

29Q d + 6 d 
2 

Rd + 2^^<i 

which yields 

af[R d a 2 ]' + af[Q d a 6 ]' = 

where 9 = 3a d /ad- These equations are valid for both averaged and non-averaged scalar quantities. 

We can see that in the case of spherical symmetry, the resulting governing equations are equations for scalar 
quantities under some circumstances. Therefore, we may be able to compare our results with the work of Buchert Q 
in this case. Certainly Buchert [6j can obtain a spatial curvature term in an FLRW background, and it is expected 
that our analysis is consistent with the work of Buchert in a deeper sense. Unfortunately, the perfect fluid case is 
very complicated for comparisons, and so it is more sensible to try to make the comparison in the dust case (i.e., the 
LTB model). 




The Dust Case 



Let us consider an inhomogeneous universe with irrotational dust. Following Buchert [6J|, one then obtains from 
Einstein's equations the following equations of motion for the effective scale factor, ap, 



3— = - T (p) v + Qv, (70) 
ax> £ 



>(^) = « 2 (P)v - \kd - \Qv I'll 



together with 



-Q t {a>vQv) +a "DQ t 



iQv) + a 4 v -(a 2 v K v ) = . (72) 



Here 1Z is the spatial scalar curvature (not necessarily isotropic) and 

Q v = l({9 2 ) v -{9)l)-{a ij a^) v . (73) 

From the analysis of the spherically symmetric models we see that, in principle, we can obtain precise evolution 
equations for TZx>, Qv (and, more generally, for the additional terms in the averaged equations for perfect fluid models). 
Indeed, since we actually take averages in our analysis we will obtain specific forms for the averaged quantities (i.e., 
in principle we will obtain explicit expressions for 1Z, Qv)- We intend to study this in the LTB models in more detail 
elsewhere [23j. 
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However, in an analysis of the spherically symmetric collapse of dust in the Newtonian regime (assuming an irro- 
tational velocity field of the form v = v(r)e r and other physical restrictions) [4l|, it was found that the backreaction, 
Qd, vanishes (as might be expected in the Newtonian approximation). Therefore, in this case the only effect of 
averaging in the Buchert approach is through a spatial curvature term. Unfortunately, the Newtonian limit is not 
expected to capture the relativistic backreaction; a GR treatment is needed to discuss the global effects of averaging. 

In our relativistic MG approach let us assume a Newtonian-like 4- velocity field of the form 



1 



[M,0,0], 



where A and B are the metric functions and v = v(r) is assumed to be small (i.e., v <C 1). 
4-acceleration is then given by 



(74) 



The corresponding 



F 



2(B~Av 2 ) 2 



B' A 



,0,0 



where 



F = AA t v 3 + (2AB r - A r B)v 2 + (AB t - 2A t B - 2ABv r )v - BB r 



(75) 



(76) 



If we assume that the fluid is pressure- free (i.e., dust), then the acceleration is zero. To lowest order in v, this 
implies that B r = 0, and from earlier we conclude that the correlation tensor is of the form of a spatial curvature. 
Consequently, in this approximation the Buchert approach (4lj | and the MG approach are consistent. 
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